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of an Expanding Bose-Einstein-Condensed Gas 
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We propose an experimental scheme to observe the Gibbons-Hawking effect in the acoustic analog 
of a l+l-dimensional de Sitter universe, produced in an expanding, cigar-shaped Bose-Einstein 
condensate. It is shown that a two-level system created at the center of the trap, an atomic quantum 
dot interacting with phonons, observes a thermal Bose distribution at the de Sitter temperature. 
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Cosmology in the early universe is a branch of physics 
which is, for all too obvious reasons, removed far from 
experiment. To use a concise formulation of the pri- 
mary dilemma which cosmology faces, there exists noth- 
ing like the possibility of "reproducing" experiments, be- 
cause there has in fact "only ever been one experiment, 
still running, and we are latecomers watching from the 
back" 0. While there can be no truly experimental cos- 
mology, as regards in particular the reproduction of the 
(presumed) extreme conditions which prevailed in the 
early stages of the universe, one might look for analogous 
phenomena in condensed matter experiments, which are 
indeed reproducible, and can be done at energy scales 
smaller by many orders of magnitude 2]. 

An archetypical model of cosmology is the de Sitter 
universe, in which space is essentially empty, and the 
curvature of space-time is due to a nonvanishing cosmo- 
logical constant A 0]. The de Sitter space-time is used 
in inflationary models of the cosmos displays a cos- 
mological horizon, and is associated with a thermal spec- 
trum of particles at a temperature Tds oc Va 0. 

Atomic Bose-Einstein condensates (BECs) 0] have 
emerged as one of the most suitable condensed mat- 
ter systems for the simulation of quantum phenom- 
ena in effective space-times II, U 1^ 1^ ■ In 
the present study, we investigate an effective H-l- 
dimensional (1-l-lD) de Sitter universe and the associ- 
ated thermal phonon spectrum of the Gibbons-Hawking 
type . We examine the propagation of axial low-energy 
modes in expanding, strongly elongated condensates, and 
find that the spectrum consists of one massless (phonon) 
mode and a sequence of massive excitations, moving in 
an effective curved 1-l-lD space-time. It is demonstrated 
that a 1-t-lD version of the de Sitter metric 0] can be real- 
ized by the massless phonon mode in a linearly expanding 
Bose-Einstein condensate with constant particle interac- 
tion. We show that an atomic quantum dot (AQD) [isj . 
placed at the center of the cloud, can be used to measure 
the Gibbons-Hawking quantum process. The detector's 
coupling to the superfluid is constructed in such a way 
that the natural time interval of the detector is equal 
to the time interval in the de Sitter metric. Therefore, 
the detector inside the expanding condensate is capable 
of measuring a thermal state at the de Sitter tempera- 



ture. The Gibbons-Hawking effect is a curved space-time 
generalization of the Unruh-Davies effect, in which a con- 
stantly accelerated detector in vacuum responds as if it 
were placed in a thermal bath with temperature propor- 
tional to its acceleration J^,il^]. Therefore, our proposal 
represents a means to confirm the observer dependence of 
the particle content of a quantum field in curved space- 
time [l^ . 

Our approach is based on the by now well established 
identity of the action of a massless scalar field propa- 
gating on a curved space-time background in D -I- 1 di- 
mensions and the action of the phase fluctuations $ in 
a moving inhomogeneous superfluid 0, 0, 0| (we set 
h = m = 1, where m is the mass of a superfluid particle): 



S = I d^'+'x— 
2.9 



-(^^-v-V^j +c2(V$)2 



(1) 



Here, v{x,t) is the superfluid background velocity, 
c{x,t) = y/gpo the velocity of sound, po{x,t) the back- 
ground density, and g = Anas is the coupling constant 
describing the short-range interaction between the parti- 
cles in the superfluid; Ug is the s-wave scattering length. 
In the second line of the conventional hydrodynamic 
action is identified with the action of a minimally cou- 
pled scalar field in an effective curved space-time. In 
the following, we derive the effective 1-l-lD action of the 
form for the lowest axial excitations of an elongated 
condensate from the full 3D dynamics of hydrodynamic 
excitations of the cigar-shaped condensate, and identify 
the metric tensor gfj,^. 

We consider a large condensate confined in a strongly 
anisotropic harmonic trap, characterized by the frequen- 
cies wii and uj±_ in the axial and radial directions, re- 



spectively {lu±_ » 



^11 J 



The initial condensate den- 



sity is given by the usual Thomas-Fermi (TF) expres- 
sion |\E'tfP = Pm(l — r'^ /R\ — z^/iJy). Here, pm is the 
maximum initial density and the squared TF radii are 
i?| = 2/i/wjj and R\ = where /x = p^g is the 

(initial) chemical potential. According to (i^l, the tem- 
poral evolution of a Bose-Einstein-condensed atom cloud 
under variation of the trapping frequencies is described 
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FIG. 1: Expansion of a cigar-shaped Bose-Einstein conden- 
sate. The stationary horizon surfaces are located at ±zh, 
respectively. The thick dark lines represent lasers creating an 
optical potential well in the center. 



by the scaling condensate wave function 

^ [r/b^, z/b) = (po + 6p)^'^ exp [i($o + *)] • (2) 

Here, $0 = - J gpoi0,t)dt + bz'^/2b + bx_r'^ /2b±_ is the 
background phase, where b^^{t) and bit) are the scaling 
parameters describing the condensate expansion in the 
radial (f) and axial (z) directions, cf. Fig.^ The mean- 
field condensate density po — \'^tf\^{'I' /b:\^, z /b) / {b\b) 
contains the scaling volume fe^fe, and 5p are the fluctua- 
tions of the density around mean-field. 

The elementary excitations above the ground state of 
the BEC in the limit LOi_/uj\\ — > 00 were first studied in 
[2ll |. The description of the modes is based on an adia- 
batic separation for the axial and longitudinal variables 
of the phase fluctuation field: 



$(r,z,t) 



(3) 



where (pnif) is the radial wavefunction characterized by 
the quantum number n (we consider only zero angular 
momentum modes). For long wavelength axial excita- 
tions, the dispersion relation reads [2l| 
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2LJl^n{n + 1) + CqA; 



(4) 



where cq = \fpj2 and k is the axial wavenumber {kR±^ <C 
1). In the action Q, we use Eq. with the rescaled 
radial wavefunction 0„ = 0n(r/6j^), and integrate over 
the radial coordinate. We then find the following effec- 
tive action for the axial modes of a given radial quantum 
number n. 



Sn = 



dtdz^ 
2ff 



+cl{d.Xnf + Mlxl] ■ (5) 



Here, the common factor Cn{z) is given by 6^C„(z) 



Ml{z,t) = gC, 



, the averaged speed of sound is c„(z, i) 



d'^rpQif)'^, and the effective mass term 
^67^ f , d'^rpnldrd'r,]'^: the radial 



cigar size is given by r„j = R\b^{l 



). The 



phonon branch of the excitations corresponds to the 
ri = solution of Eq . I^J , for which the radial wavefunc- 
tion (j)Q — const. |2l|; hence the mass term Mq ^ and 
Co(0) = 7ri?i, c(0,i) = c„=o(0,i) = co/{bA_b^/^). 

The action jSJ can be identified with the action of a 
minimally coupled scalar field in 1-1- ID, close to the cen- 
ter of the condensate. The line element reads ^3 



ds^ = Ac [-(c^ - vl)dt^ - 2v^dtdz + dz^ 



(6) 



where Ac is some arbitrary (space and time dependent) 
conformal factor and we set c = c{z = 0, t). The actions 
(0 and (O can be made consistent if we renormalize 
the phase field according to $ = Z$ and require that 
Z^b\CQ{Q)/g = l/c(0), leading to 



B 



P 



const. 



(7) 



The constant factor Z in the above relations does not 
infiuence the equation of motion, and corresponds to a 
renormalization of the field amplitude, $ = Z<^, which 
determines the coupling of l> to a detector . 

We require that the space-time metric is, in appro- 
priately chosen coordinates, identical to that of a 1-l-lD 
de Sitter universe. We first apply the transformation 
c^dt = c{t)dt, connecting the laboratory time t to the 
time variable t. Defining Vz/c — y/Hz — {Bb/co)z (note 
that the "dot" on b and other quantities always refers to 
ordinary laboratory time), this results, up to the confor- 
mal factor Ac, in the line element ds^ = —Cq{1— Az^)dt^ — 
2cozVXdidz + dz"^. We then employ a second transfor- 
mation Codr = codi + z^/Kdz/{l — Az^), with a time in- 
dependent A, to obtain the 1+lD de Sitter metric 



ds'^ 



Az^) dT^ 



1-Az 



dz^ 



(8) 



The quantity A = B^b^ /cq is a "cosmological constant," 
provided b does not depend on t. The transformation 
between t and the de Sitter time t, at a given coor- 
dinate z, then involves the exponential "acceleration," 
t/to = exp[_B5T], where ~ 27r/w|| is set by the ini- 
tial conditions. The (conformally invariant) temperature 
associated with the effective metric ^ is the Gibbons- 
Hawking temperature j5j 



ZTT ZTT 



(9) 



and the horizon(s) are located at z = ±zh = 
±i?||(w^/2^A)^/^. Combining the latter relation with 

we see that ZH/i?|| ^ 1, if Tds ^ a;||/47r. We are thus 
justified in neglecting the z dependence in Cq and c to 
describe the physics inside the horizon surfaces. 
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According to Eq. ((S)), the equation of motion 
5So/5xo = is given, at constant B, by 



1 



Co{zb) 



dz, [c^{zb)Co{zb)dz,xo] = 0, 
(10) 

where Zf, — z/bis the scahng coordinate. Apart from the 
factor Co(zf)), stemming from averaging over the perpen- 
dicular direction, this equation corresponds to the hy- 
drodynamic equation of phase fluctuations in inhomoge- 
ncous superfluids [isj . 

At t ^ — cxD, the condensate is in equilibrium and the 
quantum vacuum phase fluctuations close to the center 
of the condensate can be written in the following form 



^ = Xo = 51 

fe 



9 



4Co(0)i?||eofe 



H.C., (11) 



where a^, a[, are the annihilation and creation operators 
of a phonon. The intial quantum state of phonons, which 
we call the "adiabatic vacuum," is the ground state of the 
superfluid and is annihilated by the operators (this 
definition of vacuum and excitations corresponds to the 
"adiabatic basis" of ^(). Assuming that the expansion 
is switched on adiabatically, the density fluctuation op- 
erator, which is conjugate to xo, then reads 



4Co(0)i?||5 



x-^a.cxp 



—I 



BP 



ikzb 



-H.c(12) 







The above equation completely characterizes the n 
evolution of the condensate density fluctuations. 

The particle content of a quantum field state in curved 
space-time depends on the observer |1^]. We now show 
that the de Sitter time interval dr = dt/bB — dt/b^^^b± 
can be effectively measured by an Atomic Quantum Dot 
[l3|. Here, the AQD is an effective two-level system 
with a time-dependent level splitting (see Fig.|2J). It 
can be made in a gas of atoms possessing two hyper- 
fine ground states a and (3. The atoms in the state a 
represent the superfluid cigar, and are used to model the 
expanding universe. The AQD itself is formed by trap- 
ping atoms in the state /3 in a tightly confining optical 
potential Vopt- The interaction of atoms in the two in- 
ternal levels is described by a set of coupling parameters 
gcd = 47racd (c, d = {a,/?}), where Qcd are the s-wave 
scattering lengths characterizing short-range intra- and 
inter-species collisions; gaa = 9, CLaa = o-s- The on-site 
repulsion between the atoms f} in the dot is J7 ~ gppl^^i 
where I is the characteristic size of the ground state wave- 
function of atoms (3 localized in V^pt . In the following we 
consider the collisional blockade limit of large [/ > 0, 
where only one atom of type (3 can be trapped in the 
dot. This assumes that U is larger than all other rele- 
vant frequency scales in the dynamics of both the AQD 
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FIG. 2: Level scheme of the "Atomic Quantum Dot," which 
is embedded in the superfluid cigar, and created by an optical 
well for atoms of a hyperflne species difi'erent from that of the 
condensate. Double occupation of the dot is prevented by a 
collisional blockade mechanism. 



and the expanding superfluid. (The value of gfjp may be 
increased using Feshbach resonances 24].) As a result, 
the appropriate collective coordinate of the AQD is mod- 
eled by a pseudo-spin-1/2, with spin- up/spin-down state 
corresponding to occupation by a single/no atom in state 

a. 

Conversion of atoms between the states a and (3 is 
caused by a laser-driven transition characterized by Rabi 
frequency and detuning A. The Hamiltonian of the 
AQD interacting with the superfluid is 

ilAQD = {-A + ,g„^[po(0,i)+MO,t)]}^^^ (13) 

+^ (po(0, t)l^) exp[i($ + $o)]^T+ + H.c, 

where we use Pauli matrix notation, i.e. Uj^- — + iuy 
We introduce the wavefunction of the AQD in the form 
-0 = V/3exp[i($ -I- $o)]|/3) + V'ala). Using the scaling 
expression (0) for the condensate wave function, we find 



.dil)p 

.dlpa 



{bB[~A + p„{g^p - g)] + SV} 0^ -I- ^^Vc 



--0/3, 



(14) 



where t is the de Sitter time, and ujq — 2Q[p„il^y^'^ is 
independent of r. The perturbation potential 



6V{T)^{g^p-g)BbiT)dp{0,T) 



(15) 



provides the coupling between the AQD and the expand- 
ing superfluid (if gap — g, higher order terms in the den- 
sity fluctuations have to be taken into account in the 
Rabi term of ifT^ l. We suggest to operate the detector 
at A{t) = {gap - g)poiO,t) = {g^p - g)p„^/{b''BH^), so 
that in zeroth order in SV we obtain an effective two- 
level system with the level splitting loq, which plays the 
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role of a frequency standard of the detector. By adjust- 
ing the laser intensity contained in the Rabi frequency 17, 
one can change Wq, and thus probe the response of the 
detector for various frequencies. 

The detector response can be calculated using the two- 
point correlation function of the coupling operator SV. 
The probabilities for the detector to "click," correspond- 
ing to excitation (P+) and de-excitation (P_) of the de- 
tector, per unit of de Sitter time, are given by [13. 



dP± 1 
— — = iim — 

dr r^oc T 



T rT 



dTdT'{SV{T)SV{T'))e 



(16) 

Using the explicit form of the detector-to-field coupling 
(fT^ . we obtain 



{SVir)dVir')) oc {p.n,alY'^ 



sinh^ [TrTds (t — t' 



(17) 

This correlation function is proportional to one charac- 
terizing a thermal phonon state at a temperature Tds in 
a condensate at rest. This means that our de Sitter de- 
tector at the center of the condensate responds to the 
scaling vacuum as if it were placed in a thermal bath at 
a temperature Tds- Indeed, substituting the correlator 
(|17|l into Eq. (|16() we find that the quantities dP± j dr at 
late times are time independent, and satisfy detailed bal- 
ance conditions corresponding to thermal equilibrium at 
the de Sitter temperature given in Eq. 0: 



dP+/dT 
dP^/dr ^ 1 + ?iB ' 



(18) 



where ne = (exp[ti;o/2ds] ~ 1)^^ is the Bose distribution 
function. We stress that the AQD observes a thermal 
spectrum, even though the quantum state of the fluctu- 
ations is adiabatically connected to the initial vacuum 
state, i.e., no excitations in the adiabatic basis charac- 
terized by the dk are created (note that the time interval 
of the adiabatic vacuum is dt/Bb^ = dr/b). The latter 
process corresponds to "cosmological" quasiparticle pro- 
duction 

[HEa, which can be neglected in our situation. 
Our experimental proposal consists in the following 
steps: (a) Preparation of a large TF condensate at very 
low temperatures, (b) Introduction of an AQD at the 
center (respectively an array of AQDs to increase the 
signal), which can be used to determine the initial tem- 
perature of the cloud, (c) Linear axial expansion of the 
condensate according to (0) . The Rabi oscillations of the 
AQD can be used to monitor the thermal distribution 
nR in Since cx (p,„a3)i/2 {uj±/nf and P± cx 

|22| , the initial condensate has to be quite dense and close 
to the quasi- ID regime. The particle density throughout 
the cloud decreases like t~^. Therefore, the rate of three- 
body losses quickly decreases during expansion, and com- 
paratively long observation times are feasible. 



The AQD as a detector is suitable to measure the de 
Sitter time interval, because it couples linearly to the 
square root of the density of the superfluid gas. In prin- 
ciple, one can construct detectors coupling to different 
powers of density or superfluid velocity, and then more 
generally experimentally study the non-uniqueness of the 
particle content of various quantum states in curved 
space-time . For example, outcoupling pairs of atoms 
by photoassociation is a means to set up a detector which 
has dr/dt (X po. Finally, we note that while in this paper 
we have studied only the n = massless axial phonon 
modes, strongly elongated condensates can also be used 
to study the evolution of massive bosonic excitations. 
Together with a natural Planck scale i?pianck ~ this 
provides the opportunity to investigate, on a laboratory 
scale, the influence of finite quasiparticle mass and the 
trans-Planckian spectrum on the propagation of relativis- 
tic quantum fields in curved space-time. 
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